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Abstract: We discuss the physics of moduli (light scalar fields with Planck-suppressed 
couplings to matter) in the case of low-scale supersymmetry breaking such as gauge mediation. 
We argue that even if the mechanism of moduli stabilization is decoupled from the mechanism 
of SUSY breaking, moduli masses will generically be parametrically related to the gravitino 
mass once the cancellation of the cosmological constant is taken into account. For low- 
scale SUSY breaking, this implies that moduli fields are light, long-lived relics that will 
generically drive the universe into a matter-dominated phase, in contradiction to standard 
BBN. We discuss two scenarios for evading this problem. The first is to consider very tuned 
supergravity potentials that can make the moduli heavy enough to decay at early times and 
reheat above the temperature of BBN. Viable cosmology can be achieved in this scenario, 
which has a population of highly relativistic light gravitinos arising from decays of moduli. 
Next, we consider the more natural scenario with light moduli. The saxion field associated 
with the solution of the strong CP problem provides the most natural candidate for driving 
thermal inflation, which dilutes the moduli. We construct an explicit model for this scenario, 
when the PQ symmetry breaking scale is larger than the messenger scale. The combination 
of the constraints on relic abundance and gamma-ray flux from decay of the moduli favors a 
particular region of SUSY breaking scale, 2 x 10 3 - 10 4 TeV. For either scenario, we find that it 
is generic for low-scale SUSY to be associated with late entropy production and accompanying 
low reheating temperatures. 
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1. Introduction 

One of the most promising scenarios for how supersymmetry (SUSY) breaking is mediated 
to the Standard Model (SM) is gauge mediation (GMSB) 1 . Low-scale SUSY breaking with 
gauge mediation, compared to high-scale SUSY breaking with gravity mediation, enjoys two 
appealing features: it is automatically flavor-blind and it is a self-consistent low-energy ef- 
fective field theory yielding collider-testable predictions without relying on any knowledge of 
quantum gravity. However, it is also well-known that the gravitino, being light and stable in 
low-scale SUSY breaking scenarios, is always a cosmological embarrassment [2]. Specifically, 
for SUSY breaking scale VF e (10, 10 4 ) TeV, one has m 3/2 £ (100 meV, 100 keV). If a grav- 
itino with mass > keV was in thermal equilibrium at early times and froze out later, 
its contribution to the energy density will overdose the universe unless some late-time dilu- 
tion is present. This constraint turns into an upper bound of the reheating temperature [3]. 
For lighter gravitinos with mass 7713/2 < keV, the overclosure constraint goes away. Yet one 
still needs to worry about possible limits from cosmic structure formation. The power spec- 
trum inferred from the Ly-a forest data together with cosmic microwave background data of 

For a review, see the classic [1] . 
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WMAP requires that the gravitino density S7 3 / 2 ^ or equivalently, m 3 / 2 < 16 eV, 

if the gravitino is a thermal relic decoupled at temperatures of the order of GeV or TeV [4]. 

The situation is exacerbated as generally there could be other very weakly-coupled light 
particles present. Examples include the massless moduli that parametrize vacuum degenera- 
cies in all known superstring theories. In this paper we will collectively refer to all light scalar 
fields with Planck-scale suppressed couplings to the low-scale theory as moduli, denoted by 
T, without referring to their origins. It is well known that such moduli fields can lead to 
serious cosmological problems, known as the moduli problem, which has received a great deal 
of attention in the context of gravity mediation [5, 6, 7, 8, 9, 10, 11, 12]. When the Hubble 
rate drops below the mass of a modulus, the modulus starts oscillating coherently around 
its minimum with Planck-scale amplitude, storing an enormous amount of energy density. 
In gravity mediation, moduli decay around the same time as the gravitino, reheating the 
universe to a temperature of Tr ~ (mr/TeV) 3 / 2 keV. A successful Big Bang nucleosynthesis 
(BBN) requires the reheating temperature to be above 5 MeV and thus the moduli masses 
to be above 100 TeV. 

Though moduli exist universally in both high- and low-scale SUSY breaking scenarios, 
they turn out to be more troublesome for the case of low-scale SUSY breaking theories, where 
the problem has received relatively little attention (though see Refs. [13, 14]). The reason is 
that generically the moduli masses are of the order of the gravitino mass if they obtain their 
masses after SUSY is broken, as is often assumed. Our first goal in this paper is to clarify 
the connection between moduli masses and the scale of SUSY breaking. In fact, there are 
scenarios in which moduli stabilization is supersymmetric, e.g., by the KKLT mechanism [15]. 
Because the dynamics stabilizing these moduli is decoupled from the dynamics breaking 
SUSY, superficially it would appear that the moduli problem is easily solved in such a setting. 
The two scales can be different, and the moduli can be made heavy enough to decay at early 
times, leaving conventional cosmology unscathed. However, the empirical fact that we live 
in a world with a very small cosmological constant (relative to the scale of SUSY breaking 
and other scales of particle physics) forces a relationship among a priori unrelated scales, and 
limits our freedom to decouple the moduli and gravitino mass scales. Without tuning, the 
moduli are light and long-lived. If the moduli live longer than the age of the universe, the 
energy density stored in their coherent oscillations exceeds the critical density of the universe, 
and leads to overclosure. If the moduli are unstable on cosmological timescales, they decay to 
two photons through the coupling TF 2 , and the observed 7-ray backgrounds place a stringent 
constraint on their density [16, 17, 18]. The various cosmological difficulties associated with 
gravitinos and moduli are summarized in Figure 1. 

In this paper, we will investigate two possible methods to alleviate the cosmological 
problems of low-scale SUSY breaking. First, it is possible, via tuning the high-scale moduli 
superpotential, to give the moduli large supersymmetric masses, making them parametrically 
much heavier than the gravitino, e.g., mj- could be above 100 TeV. They would be unstable, 
decaying to gravitinos, SM particles, and superpartners at a temperature above 5 MeV, as 
required by BBN. Primordial gravitinos will be diluted while the moduli oscillate, and if 
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Figure 1: Major cosmological constraints for light gravitino (m 3 / 2 G (100 meV, 100 keV)) and 
moduli (rriT G (10 eV, 10 MeV)). For illustration, we assume that moduli is stabilized by KKLT-type 
mechanism and the moduli mass is 100 times as large as the gravitino mass. Yet it should be borne 
in mind that the exact relation between the two masses is model dependent. 

the reheating temperature is below a TeV, gravitino overclosure will not be a problem [3]. 
However, moduli decays T — > 2-0 3/ / 2 can be a new source of additional relativistic gravitinos, 
which may spoil the success of standard BBN. We will discuss this decay in some detail. 

In the case that not all moduli are tuned to be heavy by some high-scale dynamics, a new 
mechanism needs to be in place to produce additional entropy to dilute the moduli significantly 
after they start oscillating (while at the same time diluting the gravitino) . A period of late- 
time mini-inflation is frequently invoked for this purpose [9]. Thermal inflation [19] provides 
an appealing example, with a light field trapped at a fixed value by thermal corrections 
driving the inflation. It requires a scalar with a very flat potential, which is already present 
for a field in gauge mediation with a strong CP solution: the saxion! We argue that it 
is natural to connect the solution to moduli problems with the solution to the strong CP 
problem. Through mixing with the SM Higgses, the saxion would decay to SM particles 
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which thermalize at a low reheating temperature without spoiling BBN. In contrast to a 
model proposed recently [20], we consider the case where the PQ symmetry breaking scale, 
typically constrained to be > 10 9 GeV, is much higher than the messenger scale of the low 
energy SUSY breaking. After thermal inflation, the moduli density scales as Trh(T to u/Tt.i.) 3 
where the saxion reheating temperature has to take the lowest value required by BBN Tr# ~ 
10 MeV to allow for enough dilution. The temperature where thermal inflation ends is set 
by soft masses and hence fixed around the weak scale T ro ii ~ 0(100 GeV - 1 TeV). Thus 
the only free parameter left is the thermal inflation starting temperature, which gives the 
most dilution at larger values of the thermal inflaton potential height, Ttj. *>/F. On the 
other hand, smaller y/~F implies lighter moduli and longer moduli lifetimes r oc F~ 3 , which 
are favored by the 7-ray data [16, 17, 18]. Thus there is a very strong and generic tension 
between achieving enough dilution of the moduli through thermal inflation, which favors larger 
values of yF, and avoiding the constraints on decaying moduli, which favors smaller values 
of \[F (and disfavors the model of Ref. [20]). Although there is enough model-dependence 
in the moduli sector that it is difficult to make completely definitive statements, we will 
argue that successfully diluting moduli, achieving a high enough reheating temperature for 
BBN, and avoiding constraints from 7-rays selects a very particular parameter range, with 
s/F « 2 x 10 3 to 10 4 TeV, f a « 10 11 to 10 12 GeV, and moduli masses in the vicinity of the 
gravitino mass rather than heavier as in KKLT. 

The paper is organized as follows. In Sec. 2, we will first review the KKLT mechanism 
and show that even if moduli are stabilized supersymmetrically, their masses remain of order 
7713/2- Then we demonstrate that tuning (separate from the tuning that cancels the cosmo- 
logical constant!) can stabilize the moduli with large supersymmetric masses. After that, the 
cosmological properties of the heavy moduli are studied. In Sec. 3, we turn to the saxion 
solution to the light moduli problem. We present the saxion properties and the parameter 
space where the thermal inflation followed by the saxion's oscillation around its minimum and 
decays would produce enough entropy to dilute the moduli density. We discuss other issues 
and conclude in Sec. 4. In the appendix, we give a brief review of the origin and properties 
of moduli fields. 

2. The Moduli Problem and the Tuning Solution 

We begin by considering some simple models in supergravity that illustrate properties of 
moduli stabilization coupled to a SUSY breaking sector. Because moduli masses arising from 
SUSY breaking are necessarily of order F/Mp (since moduli couplings are Planck-suppressed), 
we are most interested in cases in which moduli stabilization is decoupled from SUSY breaking 
at leading order, which at least offer the possibility of lifting the moduli masses to much higher 
scales. The two models we discuss appeared in Refs. [21, 22], but there the emphasis was 
on gravity mediation and the relevance to the moduli problem was less apparent. Here we 
discuss the results with an emphasis on scenarios with very light gravitinos. 
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We set Mp = 1 in most intermediate steps but put it back to clarify results where 
numerical orders of magnitude are important. We use the reduced Planck mass Mp = 
2.4 x 10 18 GeV, for which the gravitino mass is m 3/ / 2 = ^ M and the critical density is 
p c = 3H 2 Mp. 



2.1 A Minimal Model: O'KKLT 

We begin with a simple example to show how the scales in the moduli sector and the low-scale 
SUSY breaking sector are tied up. We couple the KKLT model [15] to the simplest possible 
SUSY breaking sector, a Polonyi model with an added Kahler potential term to lift the flat 
direction. This model is well-studied under the name O'KKLT [21], and closely related work 
appears in Refs. [23, 24]. The model is 

Kq'kklt = -31og(T + Tt) + tfX - , 

Wo'kklt = W + A 3 e~ bT + fX. (2.1) 

We will assume A ~ Mp (or another very high scale). The constant Wq originates from 
fluxes in string compactification which can be tuned to make Wq <C Mp. The dimensionless 
coefficient b is 2tt in the case of D3-brane instantons and 2tt/N if it arises from gaugino 
condensation on D7 branes. For a discussion of the physics behind the logarithmic form 
of the Kahler potential, see Appendix A. The X sector is just the Polonyi model of SUSY 
breaking with a Kahler correction that stabilizes {X) at the origin. The cutoff scale M in the 
X sector is related to the underlying strong dynamics that breaks SUSY, and could be much 
smaller than Mp. 

As the two sectors are only connected by gravity, we can find the minimum of the potential 
by first inspecting the two sectors separately. The KKLT sector possesses a supersymmetric 
AdS minimum at large modulus value 1 

D T W = W T + K T W = 

-*=- w -(-^) + ^i£r (2 - 2 » 

with potential depth 

VAdS --^Mjt* Mj- (2 ' 3) 

After SUSY breaking, the X sector will contribute to the potential Vx ~ |/| 2 - The vacuum 
energies of opposite signs from the two sectors have to be balanced against each other, / ~ 
| Wo | /Mp, to obtain an almost vanishing c.c. The modulus mass is determined by Wq, and 
thus is of the same scale as the AdS curvature and the gravitino mass. More quantitatively, 



1 Here W_i denotes the branch of the Lambert W-function (implicitly defined by z = We w ) which is 
well-defined for small negative arguments > — 1/e and takes values < — 1 [25]. 
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perturbing around the minimum of two disconnected sectors, one finds that, up to corrections 
scaling as 1/ (£>**), 

V3\W\ 



V\ m in ~ => f 



M P 

™3/2 = e G/2 = —L , (2-4) 

2vW M P 



m T = e G / 2 G TTt G TT ~ m 3/2 log , (2.5) 

where the total Kahler potential G = K + log |VF| 2 . Indeed the modulus is heavier than the 
gravitino, but only by a logarithmic factor. A recent study of supersymmetric sigma models 
in AdS has shown that generically, a theory with many KKLT-like moduli with potential 
dominantly arising from one exponential in the superpotential will have all moduli masses 
proportional to the AdS curvature scale, times a logarithmic factor [26]. 

2.2 A Tunable Superpotential 

Now we will demonstrate that the modulus 
can be made parametrically much heavier than 
the gravitino through tuning the modulus su- 
perpotential. The model appears already in 
Refs. [21, 22]. Keeping the same Kahler po- 
tential Ko'kklt, consider a racetrack type 
Figure 2: An example KKLT potential V(T), superpotential with two exponential terms in 
for parameter choices leading to a minimum at the moduli sector (plus a constant): 




U = 56.8 with depth V{U) = -1.4 x KT 14 . 



Wtunable = W + A 3 e~ bT - A'V 6 ' 7 . (2.6) 



Before discussing the detailed physics, let us outline the general strategy: we want to have 
a vacuum where the moduli are much heavier than the gravitino. Moduli masses come from 
derivatives of the superpotential W at the minimum, while the gravitino mass comes from 
(W). Thus, we want to tune so that at the minimum (W) is small compared to the terms 
that make it up. One way to do this is to first tune to produce a supersymmetric Minkowski 
vacuum, then perturb it to produce a supersymmetric AdS vacuum for the modulus, then 
couple to X to uplift. 

All three terms in the superpotential are taken to be comparable, Wo, A 3 , A' 3 ~ 0{Mp). 
If any of them is much smaller than the others, it would not influence the AdS minimum, 
the modulus mass would still be of order the AdS curvature, and the physics would resemble 
that discussed in the previous subsection. The addition of the third term of the same order 
allows us to tune the parameters Wq, A, A', b, and b' to produce a supersymmetric Minkowski 
minimum by setting DjW = and W = 0, 
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«p-*-(£)*H)- 

For a controlled solution, the modulus should still be stabilized at large volume i* S> 1. 
This could be achieved, e.g., by making A relatively larger than A'. Notice that this is quite 
different from the KKLT scenario, where the large volume stabilization is obtained through 
tuning Wq to be small. Another crucial difference is that though the gravitino is still massless, 
the modulus has already obtained a large supersymmetric mass 



2ti 0) A3 



'»r- V 3^ 2 " e-^\b-b'). (2.9) 



Next we couple the racetrack sector to the SUSY breaking sector fX with / <C m^Mp. 
However, it is easier to first ignore X and build a supersymmetric AdS minimum for the 
pure modulus sector by shifting the parameter Wq so that at the new minimum 3 \ W\* = f 2 , 
while the condition D?W = is preserved. This leads to the following constraints on the 
parameters: 

.3 -bU \'3„-b't* _ f M P 



W + A d e- M * - A V 



-6A 3 e- w * + 6'A'V^ = ^IMZ. (2.10) 

2^ 

Thus, the modulus sits in a supersymmetric AdS minimum with the AdS curvature de- 
termined by the perturbation. Compared to Eq. 2.7, the position of the minimum shifts, 
between the SUSY Minkowski minimum and the SUSY AdS minimum, but only by a sup- 
pressed amount 5t* ~ f/(Mpmr) and the solution W^ shifts to 

At the AdS minimum, the gravitino mass, of order AdS curvature, is also set by the pertur- 
bation /, 



J_ 

V3M P (2i*) 3 /2 



m ;i/2 = ^d— 7 —^ i (2.12) 



which is parametrically much smaller than the modulus mass in Eq. 2.9. Having constructed 
this AdS minimum, we can then add the X field to uplift it to a Minkowski minimum, which, 
similarly to the uplifting of KKLT in the previous subsection, has only a small effect on the 
value of t* at the minimum, as the two sectors are nearly decoupled. 

Before ending this section, we would like to summarize the differences between the tuned 
case and the KKLT case in the previous section. In the tuned case, iWjMnaMel is much smaller 
than each of its individual terms. The modulus mass is determined by the supersymmetric 
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parameters, dissociated from the SUSY breaking order parameter / that sets the gravitino 
mass. From Eq. 2.11, one could see that the tiny c.c. is achieved by a fine tuning to the order 
0(fMp/A s ) among the large Planck-scale parameters in the moduli sector. The modulus 
potential is shallow but steep as shown in the right panel of Figure 3. We stress that the 
tuning that makes the modulus heavy is different from the tuning that cancels the c.c; even 
fixing the depth of the supersymmetric AdS minimum that we uplift, the typical potential 
will look like the left-hand plot in Fig. 3, and the modulus mass will be of order m 3 /2- 




Figure 3: Two example potentials V(T), untuned at left and tuned at right, arising from the super- 
potential 2.6. As in Fig. 2, there is a minimum at t* = 56.8 with depth V(t*) = —1.4 x 10 -14 in each 
case. In the plot at right, the modulus at this minimum is tuned to be very heavy. 

Our argument that moduli masses will be of order the gravitino mass unless there is a 
tuning such that the individual terms in W, at the minimum, nearly cancel is very similar to 
the one already given in Ref. [7]. However, they explicitly state an assumption that moduli 
fields are massless in the absence of SUSY breaking. We emphasize the key point that we need 
not assume that the moduli masses arise from SUSY breaking effects to make this argument. 

2.3 Moduli Decays 

If we tune to have a light gravitino and heavy moduli, one potential problem is an over- 
abundance of gravitinos arising from moduli decays. The heavy moduli could either cascade 
through SM superpartners ending in gravitinos and SM particles, e.g., T — > 2g —)•••• — >• 
2^3/2 + • • • or decay directly to a pair of gravitinos T — > 2^3/2- The decay widths scale as 

= amy/(47rAfp) with a model-dependent constant a. Below we will calculate the partial 
widths and the branching fractions of different decay channels. 

The simplest way to couple a single modulus to the MSSM is through the gauge kinetic 
function, J d 2 9TW a W a , which leads to, in components, 

C = -—±= l^G%G a ^ + ^G%G a ^ - ^(TXX + h.c.)) , (2.13) 

where the subscripts R(I) denote the real (imaginary) component of T, and T = \J G TT \ (T — 
t*) is the canonically normalized modulus field. For K = — 31og(T + T'), a straightforward 
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calculation yields the partial widths 



AT 777^ 

r(T R>I -> gg) = T(T RJ -+ XX) = — ^, (2.14) 

where A?" counts the degrees of freedom of the final states. For instance, N = 8 for gluons 
and gluinos. The decay widths to the gauge bosons and gauginos are equal as SUSY breaking 
effects are suppressed by m 3 / 2 /w-T or m so /(/mr. Notice that the result is independent of 
the gauge coupling or the modulus VEV, which is generically not true when the MSSM is 
coupled to multiple moduli. Other couplings such as J d i 9Q^e~ 9 ^ v Q would induce three- 
body decays such as T — > qqg but these would be suppressed by the gauge coupling as well 
as a three-body phase space factor. Thus we will neglect them for the rest of the discussion. 

The direct decay of moduli to gravitinos, T — > 2ip 3 / 2 , is one potential source of cos- 
mological difficulties [27]. The coupling here depends on the SUSY-breaking .F-term of the 
modulus field. The generic expectation is that for heavy moduli Gf ~ m^n/mT- Let us 
briefly review the argument. Suppose, around the supersymmetric vacuum, we had a large 
mass, \mf{T — i*) 2 . After coupling to the SUSY-breaking field, we have a new term in the 

potential, e K K x x |/| 2 . Because e K is T-dependent, the VEV of T will shift, between the su- 

I f\ 2 

persymmetric AdS vacuum and SUSY-breaking Minkowski vacuum, by 5t* ~ This leads 

to a new F-term for T, F f ~ St* d f {D f W) ~ £k Thus, we conclude that Gf ~ ^3/2 / m r, 
in which case the decay width is order mf/Mp [27]. 

However, an important subtlety arises in calculating this decay [28] (see also the general- 
izations in Ref. [29]). The physical Goldstino is a linear combination of ipx an d ipT while the 
orthogonal combination ^ is massive. Similarly, its massive scalar partner $> is also a linear 
combination of X and T. It was first pointed out in [28] that in certain classes of models, such 
as KKLT, the mixing is supersymmetric at 0{m^jijmT). Consequently the decay amplitude 
of the modulus, or more exactly of the massive scalar $, starts at 0((m 3 / 2 / m T) 2 )- In the 
KKLT model, the two small numbers 777,3/2 /tot and 1/t* are of the same order, but this not 
true in more general models like that of Sec. 2.2, for which 1713/2 / m T C 1/t*. In order to 
disentangle which small factor appears in the general form of the argument of Ref. [28] , and 
for completeness of the discussion, we will go into more details of the argument below. 

The coupling of a heavy scalar $ to the gravitino is e G / 2 G$§'ip fl a flu ^ u /2. Thus, to 
calculate the width of one needs to work out In terms of the G function, the potential is 
V = e \G GiGtf — 3), where i runs over all the fields. To have vanishing c.c, G %% GiG^ —3 = 
0. The calculations below will assume a factorizable Kahler potential K = — 31og(T + T') + 
/c(|X| 2 ) where A:(|X| 2 ) is a function of \X\ 2 . From the stationary condition dV/dT^ = 0, one 
obtains 



G XX 'G 



Gt = drG = - ^" xtrt Gx (1 + 0(m 3/2 /m T )) 



V2 G xxt G x trtG x (1 + 0(m 3/2 /m T )) , (2.15) 



rriT 
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where we employ my = 

= m Z j 2 G TT ^ Gtt in the second line. Next we will 
diagonalize the scalar mass matrix. The off-diagonal term would induce a mixing between T 
and X, 



^§^/§S^, (2.16, 
where X is canonically normalized. Thus the heavy mass eigenstate is, assuming my > mx, 

$ = eX + T, (2.17) 

where 



m\ - m\ G Ttrt V G TTi 



m T m z/^ G xxi G TTi GxW . (2.I8) 



The factor relevant for the decay, G$, is (at leading order in m^^/mx) 
G$ = eG^ + G f = eG x + VG TTt G T 

= ( -^-^ - l) ^V^G^G^Gx (2.19) 



= S^( f x 2 )^G XX \ (2.20) 
\ mi, — m x J mx 

In the limit mx <C tot, there is a suppression of the direct decays. The argument of Ref. [28] 
has a "chirality-suppressed" coupling G<j> ~ {1^/2/ m T) 2 in KKLT coupled to a flat pseu- 
domodulus X, matching the mx —> limit of Eqn. 2.20. However, in the low-scale SUSY 
breaking limit we are interested in, where the mass of X is present in effective field theory in 
the Mp — > 00 limit, we will generally have m x /m^ 3> m^^jm^ and Eqn. 2.20 is the domi- 
nant term. For generic low-scale SUSY breaking, with completely general Kahler potential, 
any suppression beyond G$ ~ m 3 / 2 / m T is absent, leading to a decay width 

3 

r(T^2V 3 / 2 ) = — (2.21) 
v vs/ij 2887TMJ,' v ' 

where c is a model-dependent order-one factor. For the nongeneric form K = — 31og(T + 
Tt) + k(\X\ 2 ), c = 27m x /(k xx ,(m 2 T - m x f). 

2.4 Cosmological constraints 

Heavy moduli have a relatively high reheating temperature 

Tr = (^ 2 5 «/90)- 1/4 v / WV 

-3^rf m T \ 3/2 



5.5 x 10~ 3 MeV (^y) , (2.22) 
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where we approximated the modulus decay width by Tt = m^/(A-KMp). For mr > 100 
TeV, Tr > 5 MeV, compatible with the requirement for a successful BBN. 2 In the range 
y/~F G (10 3 — 10 4 ) TeV, m 3 / 2 G (1 — 100) keV, the decay of the heavy moduli with mass rriT G 
(100 — 10 6 ) TeV could reduce the gravitino density below the critical density. However, one 
still needs to worry about the gravitinos from the direct decay of the moduli T — > 2ip3/2- The 
energy they carry could increase the Hubble expansion rate, resulting in an overproduction of 
He 4 . More quantitatively, the energy density of the gravitinos has to be smaller than 20% of 
that of the SM particles, or equivalently, Br(X — > 2ip 3 / 2 ) < 0.2. Combining Eq. 2.14, 2.21, 
for the tuned model we consider, this constraint turns into c < 18, which could be easily 
satisfied for mx < in the case of a nongeneric Kahler potential. For a completely generic 
Kahler potential, however, this moduli-induced gravitino problem is a dangerous constraint 
on the tuned scenario [27]. 

One might also worry about whether similar problems would be caused by gravitinos 
at the end of the cascade decay chains of moduli. In this case, however, the superpartners 
produced from the decaying moduli would first enter the thermal bath, leaving no chance 
for them to decay to gravitinos as the thermalization time scale C(GeV _1 ) is much shorter 
than the lifetime of the decay to gravitino ~ F 2 /m 5 so j t 3> 0(GeV~ 1 ). Yet, eventually the 
next-to-the lightest superparticles (NLSPs) will freeze out, and their late decays to the grav- 
itino after the BBN epoch would modify the abundances of the light elements and thus are 
constrained [3]. For instance, for f as NLSP, its lifetime has to be T _1 < 6 x 10 3 sec [31]. 
This is easily satisfied for low-scale gauge mediation. 

We close this section by noting that, if all moduli are tuned to be heavy, any QCD 
axion in the theory should not be a fundamental (e.g., string theory) modulus, but rather a 
conventional field coupling with larger than gravitational strength. In particular, because the 
QCD axion should get its potential dominantly from QCD instantons, in the limit of unbroken 
SUSY the QCD saxion will also be a flat direction, and the saxion mass scales as F/Mp if it 
is a fundamental modulus. Thus, the QCD saxion can never be tuned to be heavy without 
ruining the solution of the strong CP problem. (General comments on axions as unstabilized 
fundamental moduli may be found in Ref. [32].) A field with couplings that are not Planck- 
suppressed, however, can play the role of the QCD axion, as the stronger couplings can allow 
the saxion to get a larger SUSY-breaking mass. Thorough recent discussions of axion and 
saxion couplings and cosmology in the limit when moduli are assumed to be decoupled (as 
in this tuned scenario) appear in Ref. [33]. However, given that any theory with low-scale 
SUSY breaking and a QCD axion should contain a saxion field that is not a Planck-coupled 
modulus, we find it more compelling to avoid tuning, keep the moduli light, and use the 
saxion dynamics to dilute their abundance. We now turn to this scenario in the next section. 

2 A related discussion of the need for tuning to make moduli heavy appeared recently in Ref. [30], which 
argued that moduli should be made heavier than the messenger scale so that they can be integrated out. We 
will not impose such strict demands; because moduli couple very weakly, they don't affect the dynamics of the 
SUSY-breaking sector in detail, so we only consider phenomenological and cosmological constraints on their 
masses. These considerations only require that they be heavier than 100 TeV. 
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3. The Saxion Scenario 



3.1 The Model and the PQ Scale / SUSY Scale Relation 

In this section we will discuss what we consider to be a natural solution to the moduli problem 
in the gauge mediation context: namely, diluting the moduli through cosmological dynamics 
driven by the saxion field associated with the solution of the strong CP problem. This will be 
a concrete example of the general paradigm of thermal inflation, as developed by Lyth and 
Stewart [19], and originally applied to low-scale gauge mediation in Refs. [13, 14]. By tying 
it to other physics, we aim to make the general picture more concrete and palatable. Closely 
related recent work by Choi et al. appeared while this paper was in preparation [20], although 
we will mostly focus on a different part of the parameter space for the model. 1 Another similar 
model, with an emphasis on the n/Bfi problem, appeared recently in Ref. [35]. Related work 
in the gravity mediation context can be found in Ref. [36]. 

We consider a model with two singlets S and Y and vectorlike quarks Qi,Q 2 with su- 
perpotential 

S 2 ( S \ n 

W = \ Q SQ l Q 2 + \H 1 ^H u H d + X Y { — \ S 2 Y + W G msb, (3.1) 

where we will focus on cases with n = 1,2. Wqmsb includes the MSSM superpotential 
(except for the /i-term, which is supplied by Ah above), interactions with messengers, the 
SUSY-breaking sector, and so on. This model is a version of the Kim-Nilles approach to 
the /i problem, linking the PQ and TeV scales [37]. In the supersymmetric limit, the |-Fy| 2 
term in the potential causes the dynamics to favor (S) = 0. On the other hand, after SUSY 
breaking, gauge mediation will affect the Qi fields, which in turn affect S. 

One desideratum of an axion model is to have 10 9 GeV £ f a ^ 10 12 GeV, where the 
lower bound comes from stellar cooling and the upper bound from requiring that the axion 
relic abundance does not overdose the universe. (See Ref. [38] for a recent review of axion 
physics.) In our model, we should make the identification: 

f a = V2(S), (3.2) 

because after integrating out the heavy quarks Qi,Q 2 , we find an induced coupling: 

[ d 2 0— log S trW a W a + h.c, (3.3) 
J oir 



and using S = ((S) + ~^) e , we find that this includes: 



tt.s 



Saxion D ^raF« u F a ^, (3.4) 

with f a = v2 (5), which matches the standard normalization of the axion-gluon-gluon cou- 
pling. 



1 A preliminary version of our work was presented earlier [34]. 
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As with any field theory model that aims to solve strong CP, the expectation that grav- 
ity breaks continuous global symmetries will require that Planck-suppressed PQ-violating 
operators added to the superpotential 3.1 or the Kahler potential are forbidden by a dis- 
crete symmetry [39, 40]. More specifically, take an additional PQ breaking operator in 
the superpotential S m+3 /M^ as an example. Its contribution to the axion potential is 
SV ~ f™ +2 /Mp\F s \ ~ f™ +3 v 2 /M™ +1 with v = 246 GeV. For f a = 10 9 (10 12 ) GeV, this 
would overcome the QCD contribution 9m 2 f 2 unless m > 3 (6), corresponding to N > 6 (9) 
for an underlying discrete symmetry. 

For concreteness, we will discuss the case of minimal gauge mediation; i.e., we will assume 
a field X, (X) = M mess + 9 2 F, with superpotential coupling XQ& to messengers in the 5 
and 5 of SU(5). Extending the discussion to the more general case [41] is straightforward. 
At the messenger scale, the soft mass m 2 s = (up to threshold corrections), whereas <5i,2 get 
positive squared masses from GMSB. In the minimal case we have (neglecting effects of RG 
running, and assuming QCD couplings dominate): 

The usual one-loop RGE for the soft masses, 

16 ^ 2 ^g7? = 3 |Aq|2 H + m Qi + m ^ + \ Ax * I 2 ) ' (3 - 6) 

shows that for (S) < M mess , RG flow drives S tachyonic and leads to a nonzero VEV that will 
be determined by an interplay between the SUSY-breaking potential (which wants S to be 
large) and the tree-level | Fy \ 2 potential (which wants S to be small) . There are two different 
scenarios for the parametric behavior of the potential, depending on the relationship between 
PQ and messenger scales: 

• The PQ scale is smaller, (S) ^ M mess . We assume M mess is not much larger than the 
PQ-breaking scale. In this case we can use Eqs. 3.5 and 3.6 to estimate 

mS[S) ~ ^^ l0g ^P (3 ' 7) 

and minimizing the potential V(S) w m 2 s S^S + |Fy| 2 gives: 



2 



(\s\) 



2n+2 _ a s 



8(n + 2)7T 4 



Ay 



/ FMfi \ 2 M mess 



log^=. (3.8) 



Mmess J Aq (S) 



Crudely, then, the relationship between the PQ-breaking scale and the SUSY-breaking 
scale in this model is: 

/ FMu \ ^ i 
V Mmess J 
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For n = 1, the Kim-Nilles \x term is automatically of the same order as the soft 
masses. As a numerical example, if we fix n = 1, Xq = Ay = 1 and th/q. = (1 TeV) 2 , we 
find that M mess = 10 12 GeV corresponds to (S) = 5.2 x 10 10 GeV. This particular choice 
corresponds to \f~F = 2.8 x 10 8 GeV and 1713/2 = 4.4 MeV. In general, this scenario favors 
777,3/2 ~ 0(100 keV) or somewhat larger, corresponding to relatively high-scale gauge 
mediation and stable NLSPs on collider time scales. Because the scenario is discussed 
more extensively in reference [20], we will focus more on the following scenario, which 
allows lighter gravitino masses. 

The PQ scale is larger, (S) ~3> M mess (for calculability, we assume a hierarchy). In this 
case the correct calculation of the potential appears in Section 5.2 of [42] as an example 
of analytic continuation into superspace. Assuming that QCD effects are dominant, 
and keeping only the largest power of log — > their result is: 

lq] dV eS _ a 2 s 2 2 \ Q \S\ 

|5| a^""i2^ |F| log m„ ' (3 - 10) 



L mess 



Minimizing the combination of this and the \Fy\ tree-level potential, we find: 



Ay M, 



2 



(3.11) 



24(n + 2)tH 

Thus, schematically, in this example we have a scaling: 

(3.12) 

so that even though the establishment of the PQ-breaking VEV depends on SUSY- 
breaking dynamics, the scale of PQ breaking is much larger than the scale of SUSY 
breaking. This allows low-scale SUSY breaking, with NLSP decays to gravitino on 
collider timescales, to be consistent with axions in the allowed window for conventional 
axion cosmology. 

We will now turn to a more detailed study of the latter case. In Figure 4, we illustrate 
the relationship between the SUSY- breaking scale, \/~F, and the axion decay constant f a . 
For n = 1, it is not possible to extrapolate to larger y/F and f a . The reason is that, if we 
want soft masses to be in the TeV regime, M mess must be taken proportional to F/m so ft, 
and for \/~F <; 10 4 TeV, M mess becomes larger than (S) and the model enters the higher 
scale SUSY-breaking regime studied by [20]. From the left panel in Figure 4, the accessible 
range of f a for n = 1 is roughly 10 9 GeV ^ f a ^ 10 10 GeV, and achieving such values 
of f a for low-scale SUSY breaking is easier at small Ay. (Of course, in general, one could 
interpolate between the two regimes; there is no physical bound on yF or f a , just a region 
in which neither approximation we have described is quite accurate.) For n = 2, the scaling 
fa ~ (-FMp) 1 / 4 allows for a larger separation between the PQ scale and the SUSY breaking 
scale. As shown in the right panel in Figure 4, our model with n = 2 probes the region 
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Figure 4: The axion decay constant f a versus the SUSY-breaking scale vF, based on Eq. 3.11. To 
fix M mess , we require that the squark soft masses (Eq. 3.5) are 1 TeV. The coupling Xq = 1 for all 
curves (because it appears only in the log, varying it has little effect.) The (from bottom to top) 
long dashed grey line, dotted orange line, solid black line, dashed purple line, dotted red line, and 
dot-dashed blue line are for Ay equal to 100, 10, 1, 0.1, 0.01, and 0.001, respectively. The horizontal 
green dashes at 10 9 GcV signal the lower bound on f a consistent with stellar cooling, while those at 
10 12 GeV mark the bound at which axions overdose the universe. 



10 10 GeV £ f a & 10 12 GeV, at the upper end of which the axion could be cold dark matter. 
In this case, large Ay ^0(1) are preferred as for fixed vF, small Ay could push f a above 
the upper limit. 

3.2 Saxion Mass and Decays 

In the low-scale SUSY breaking regime discussed above, it is straightforward to calculate that 
the saxion mass around its minimum is: 



12vr 4 (S) 



, (n + 2) log 2 



Aq (S) , X Q (S) 



Mn 



log 



Mr, 



(3.13) 



As illustrated in Figure 5, the resulting masses can range from MeV to 100 GeV for the model 
choices that give us low-scale SUSY breaking with reasonable axion decay constants. Because 
couplings are 1/ f a suppressed, there is no bound on such particles even when they are light. 
Saxions always have a decay to a pair of axions, which arises from the kinetic term: 



r( s 



m" 



aa) 



327T/2 



(3.14) 



This will lead to a disaster for cosmology if it is the dominant decay mode: because the saxion 
dominates the energy density as it oscillates coherently around its minimum, if it decays 
predominantly to axions the universe will thereafter be full of highly relativistic particles 
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Figure 5: Saxion mass versus the SUSY-breaking scale VF, based on Eq. 3.13. Parameters are as in 
Figure 4, with (from top to bottom) long dashed grey line, dotted orange line, solid black line, dashed 
purple line, dotted red line, and dot-dashed blue line are for Ay equal to 100, 10, 1, 0.1, 0.01, and 
0.001. Points are plotted only if 10 9 GeV < f a < 10 12 GeV for the given parameter choices. 



that do not thermalize with SM degrees of freedom. We must require that the dominant 
decay of the saxion is to SM fields. Note that the decay to two gluons from the operator 
SW a W a does not help, as the associated partial width is smaller by a factor of relative 
to the partial width to two axions. 

This is where the S 2 H u Hd superpotential term plays its key role, providing the crucial 
mechanism for saxion decays to lead to reheating. The mass mixing term for Higgs and saxion 
(taking the decoupling limit for the MSSM Higgses) is: 

C D -j^ysh. (3.15) 

(In particular, tan/3 dependence drops out.) For instance, if m s > 2rrib, we have a decay 
s — > bb with partial width: 

For smaller m s , e.g., 2m c < m s < 2m&, the saxion could decay to charmed mesons (or r's). 
For n ^ 2rrih, the saxion would decay to the SM 99% of the time and thus the energy density 
contained in relativistic axions does not spoil BBN. Near the upper bound of f a , 10 11 GeV 
^ fa & 10 12 GeV, the region that our model with n = 2 covers, the reheating temperature, 
Trh-, is lower, ranging from a few MeV to 100 MeV. On the other hand, near the lower bound 
of f a , 10 9 GeV ^ f a ^ 10 10 GeV, the region that our model with n = 1 covers, Trh could be 
as large as O(10) GeV. 
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3.3 Cosmology 

We wish to dilute the abundance of moduli fields, which we assume have a mass my of 
approximately the order of m 3 / 2 (but perhaps somewhat larger). Specifically, we will write: 

mr^^ p , (3.17) 



™3/2 

matter-dominated for 3H ^ m^. As explained in Ref. [9], at this point we have: 



where for instance in KKLT ^ ~ log — — can be ~ 10 2 . The universe rapidly becomes 



T - 2 "^° -0.05-9=, (3-18) 



where we take the effective number of degrees of freedom to be g* s = 275 and the moduli 
oscillation amplitude to ~ Mp. rij,/s ra _d remains fixed in the absence of entropy production. 
The temperature of the radiation at the point when the moduli dominate the energy density 
is T m .d. ~ (90/(vr 2 ^ s )) a25 v / P" » OAy/^F. Due to the t /y/F factor, for low-scale SUSY 
breaking y/~F ~ (10 — 10 4 ) TeV, n^/s can be as large as (10 9 — 10 12 ). Of course, the number 
that must be small at late times is Pf/s, which shrinks with my, and can be estimated to be 
~ (10 3 — 10 6 ) GeV. Given the critical density of the Universe versus the present entropy is 
p c /s = 3.6 x 10~ 9 /i 2 GeV, where h is the Hubble rate, the moduli density has to be diluted 
by a factor of ~ (10 12 - 10 15 )! 

Let's consider how the saxion can help to solve this problem by realizing the idea of 
thermal inflation [19]. By integrating Eq. 3.10, and assuming the cosmological constant is 
approximately zero at the PQ-breaking zero-temperature vacuum, we see that the vacuum 
energy when S = is approximately given by: 

V " S ViS = °> " 3^ |F|2 ** £t. ~ ( (M - " A) • (3J9) 

where the range of numbers quoted is extracted from scanning the set of model lines displayed 
in Fig. 4. Thus, we expect the vacuum energy for S at the origin corresponds to a scale 
slightly, but not far, below \ ~F. Thus, when the Hubble scale drops to order y/F, this constant 
energy density provides a potential source of late inflation. Because we have U 1/4 ~ 0.1 VF, 



and Tm.d. ~ V F, moduli domination and the onset of thermal inflation happen almost at 
the same time, but by assuming that £ is somewhat large we will simplify the discussion and 
speak as if events are well-separated in time. Thermal inflation begins when Vq dominates 
the energy density, i.e. Vq > m^rij,, which occurs when the radiation is at a temperature: 

30U x 1/3 



T T .i. w „ " ~ (0.02 - 0.06) VF. (3.20) 

V 71 " 9*s-Lm.dJ 

In order for Vo to lead to inflation, we need S to remain at the origin for a period of 
time. But this happens automatically in our model: at high temperatures, the saxion field 
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will be trapped at the origin by thermal corrections to its potential. These corrections arise 
from interactions with the Q fields, which keep it in thermal equilibrium with SM degrees of 
freedom. In particular, the leading thermal term in the potential goes as [43]: 

^hermal(5) = ^|Aq| 2 T 2 |5| 2 . (3.21) 

Thus S remains at the origin until the universe cools enough that the tachyonic mass from 
eq. 3.6 becomes larger. Then, S begins to roll out to its VEV when the temperature drops 
to a scale set by the soft masses of the Qs: 



-roll 




0.5 to 1 TeV. (3.22) 



During the thermal inflation, the number density of f redshifts as e 3Ht while the tem- 
perature T redshifts as e~ Ht . Thus we have (using temperature as a clock): 

T TO \\ \ 3 



n f (T roll ) = n f (T TJ .) ( ^ j . (3.23) 

In particular, the ratio of pf = mrrif to the total energy density changes by this factor, 
because the total energy density remains approximately Vq throughout thermal inflation. 

After thermal inflation, the energy Vq is converted to energy mgng for the saxion s 
oscillating around its minimum. The ratio rif/ns remains fixed as s oscillates: 



m s Pf m s I iron 



n$ rriT ps m T V^T.I. 
When s decays, the entropy changes to 



(3.24) 



_ 271-2 ^3 _ 4 Prad _ 4m 5 n 5 
Srad - ^9*sT RH - -— - - — , (3.25) 

with Trh the reheating temperature of s. From this one can conclude that the final value of 
the energy density in the modulus versus the entropy, in terms of the temperatures, is: 

Pf nf n s m s ( T ro n \ 3 ST RH 3 { T Ion \ 3 

= m T m T _ — ^ -T RH - — , 

Srad n s Srad m T \T T jJ 4m s 4 \Ttj.J 

2 Pf/ Sr&d n -, n 8 Pf _ w n 8 T RH ( ^roll 



, = rrip^t « 2 x 10 s « 2 x 10 s ^ [=^ L \ (3.26) 

Pc/Srad Srad CxeV \T T jJ 

Thus, we want to lower Trh while making the ratio T r oii/^T.i. as small as possible. It turns out 
that for relatively large SUSY breaking scale y/~F ^ 2x 10 3 TeV and low reheating temperature 
Trh ~ 10 MeV, the saxion oscillation and decay could reduce the moduli density to be below 
the critical density. This implies that only our model with n = 2 works, where f a could be 
close to 10 12 GeV, hence the saxion mass is relatively small, and consequently Trh could be 
as low as 10 MeV. 
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As briefly mentioned in the introduction, for a modulus with mass in the range 200 keV 
to 10 GeV, even if its lifetime is longer than the age of the Universe, its decays to photons in 
the past would contribute to the 7-ray flux [16, 17, 18]. Thus the observed 7-ray background 
provides a much more stringent constraint on the moduli density, as shown in Fig. 1 in [17]. 
The bound gets strongest for moduli with mass around 100 MeV, requiring n f h 2 < 10~ 10 ! 
As this constraint is obtained by assuming that moduli dominantly decay to two photons 
through the operator cTF^F^ with c a model-dependent coupling, it might be relaxed if c 
is suppressed, e.g., by the gauge coupling c ~ a em j (47r), because the moduli lifetime becomes 
longer and fewer have decayed at a given time. Considering that the ultraviolet theory could 
possess many moduli, one would still need to worry that at least one modulus will have 
unsuppressed width to photons. With this in mind, assuming this constraint applies (at least 
to some of the moduli), one would find several interesting implications. First, in the KKLT 
moduli stabilization, the modulus mass is fixed at roughly 100 times the gravitino mass; for 
\[F ^ 2 x 10 3 TeV, mT J> 200 keV, and thermal inflation, though capable of reducing the 
moduli density to just below the critical density, could not dilute it sufficiently to satisfy the 
constraint from the 7-ray flux. One has to invoke a more efficient dilution mechanism, or 
adopt a model of moduli stabilization with a smaller ratio of modulus to gravitino mass. This 
constraint also applies to relatively high scale mediation with yF > 10 4 TeV as considered 
in [20]. The region of parameter space under the least tension from the various constraints 
has y/~F ~ 2 x 10 3 to 10 4 TeV, Trh ~ 10 MeV, and moduli masses ~ ^3/2 without the 
log-enhancement as in KKLT. 

4. Discussion 

We have argued that, generically, low-scale SUSY breaking suggests the existence of light 
moduli fields with masses approximately of order m 3 / 2 . One of the most interesting impli- 
cations is that the universe likely had a period of late entropy production. If we evade the 
generic scenario by tuning the moduli to be heavy enough that their decays reheat the uni- 
verse to above 10 MeV and allow for BBN, the entropy production arises from the decay of 
the moduli fields themselves. If we aim for a more natural scenario, the entropy production 
must come from some other source; we have suggested the saxion as a natural candidate. In 
the case of gravity mediation with to 3 / 2 of weak scale, which we have not considered here, 
one can similarly argue that a reasonable scenario is that moduli are at around 100 TeV and 
their decays reheat the universe (for several incarnations of this idea, see Ref. [44]). 

The big picture is that considerations of moduli cosmology generically lead us to expect 
an association between SUSY and low reheating temperatures, for any scale of SUSY-breaking 
that could lead to weak-scale superpartners. Another line of reasoning, completely indepen- 
dent of moduli problems, has recently led to a similar conclusion. Namely, the combination 
of gravitino cosmology [3] and the cosmology of a QCD axino, considered in tandem, require 
low reheating temperatures [45]. 1 

1 This follows earlier work on axino cosmology, including the possibility of axino dark matter and bounds on 
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One exciting possibility suggested by these considerations is that there is a great deal 
of physics happening at scales that are not far away. This potentially extends beyond the 
saxion physics and the SUSY-breaking sector. For instance, given low reheating temperatures, 
one might expect that baryogenesis happens at low scales, even below the electroweak phase 
transition, and perhaps can be probed at colliders. Scenarios such as those of Ref. [51] 
potentially fit into the larger picture of cosmology that we have sketched. An alternative 
possibility, advocated in [13], is Affleck-Dine baryogenesis that produces such a large baryon 
asymmetry that it can explain the observed number even after being diluted by late entropy 
production. However, such a scenario deserves a careful check in light of possible Q-ball 
dangers [52]. 

Though we have focused mostly on the cosmology of low-scale SUSY breaking, we want to 
mention a possible collider signature related to the saxion scenario. There in some cases (e.g., 
n = 1 models), the superpartner of the saxion, the axino, has mass at or below the weak scale, 
rria ^ and mixes with the Higgsinos through the superpotential operator S 2 H u Hd. Thus 
it's possible to have the NLSP (e.g., Higgsinos) decaying to axino plus Z or Higgs with collider 
length scale [53, 54]. The signature can mimic signatures previously studied in the GMSB 
context [55, 56, 57]. However, for the axino decay mode to be comparable to that to gravitino, 
the PQ scale has to be relatively low, f a ^ 10 10 GeV, which doesn't work well for the saxion 
being a late-time entropy production source. Yet this collider signature could well be present 
for relatively high scale SUSY breaking and mediation. It raises the interesting question: 
suppose we find evidence for a neutralino that propagates over a macroscopic distance and 
decays to a Z or a Higgs boson plus an invisible light particle. Should we interpret it as a 
decay to a gravitino, an axino, or something else? Such decays offer the prospect of measuring 
important scales like y/F or f a , but only if the correct interpretation is found. 

In summary, we see many interesting possibilities for consilience between TeV-scale par- 
ticle physics, including accessible collider signals, and the fundamental physics of the Planck 
scale, driven by the cosmological problems posed by very weakly interacting particles. It 
would be interesting to pursue more complete models, including cases with multiple moduli, 
and to compute the cosmology of saxion-driven thermal inflation in more detail. 

Acknowledgments 

We thank Ian- Woo Kim, Masahito Yamazaki, and especially Josh Ruderman for involvement 
in an early stage of this work. We thank David Shih for useful comments on the manuscript. 
We would also like to thank Nima Arkani-Hamed, Zackaria Chacko, Nathan Seiberg, and 
Tomer Volansky for useful comments, suggestions, and criticisms. J.F. is supported by the 
DOE grant DE-FG02-91ER40671. M.R. thanks the PCTS for its support. L.T.W. is sup- 
reheating temperature [46, 47, 48]. See Ref. [49], however, for an argument that in particular models, especially 
of the DFSZ type, this bound might be circumvented. Further recent work on cosmology of gravitinos, axinos, 
and saxions appears in Ref. [50]. 



- 20 - 



ported by the NSF under grant PHY-0756966 and the DOE Early Career Award under grant 
DE-SC0003930. 



A. Review: Origin and Properties of Moduli Fields 

In string or M theory, one finds that parameters like gauge couplings are never fixed numbers, 
but are determined by the expectation value of scalar fields. For instance, if a gauge field in 
the 4D effective theory at low energies arises from D7 branes wrapping a 4-dimensional cycle 
S in a Calabi-Yau manifold, the low-energy gauge coupling is determined by the volume 
of that cycle, which is a dynamical quantity because the theory includes gravity. More 
precisely, \ oc Vol(S). In a supersymmetric theory, this is the real part of a chiral superfield 
with an imaginary part realted to the 6 angle of the theory, 6 oc J* 2 C4, where C4 is the 
Ramond-Ramond 4-form of type IIB string theory. This is a representative example of a wide 
class of moduli in string theory which have approximate continuous axionic shift symmetries, 
9^0 + 58, broken to discrete symmetries (e.g., 6 — > 9 + 2ir) by instantons. Examples are 
Kahler moduli in IIB strings and both Kahler and complex structure moduli in IIA strings. 
Moduli that determine the values of gauge couplings are always of this form, because there 
is always an associated #-angle. This means, in particular, that any superpotential terms 
involving such moduli are of the form Ae~ bT , with possible values of the number b constrained 
by the remaining discrete shift symmetry. 

Because the moduli arise from closed string modes like the graviton, Ramond-Ramond 
fields, or B field, they couple to normal matter with gravitational strength (i.e., through 
1/Mpz-suppressed operators). While we know that our universe doesn't contain massless 
scalars with gravitational-strength interactions, fifth-force searches still allow fairly light fields 
with such couplings to exist. At the very lowest possible scales of SUSY breaking, yF ~ 10 
TeV, direct searches for these forces may be relevant [58], with current experimental bounds 
excluding gravitational-strength forces at a range A = 56 //m = (3.5 meV) -1 [59]. 

One other property of moduli that will be of interest for us is that they generally have 
noncanonical Kahler potentials that diverge as moduli approach certain special points. For 
instance, in the D7-brane gauge coupling example, the limit where Vol(S) — > is associated 
with the existence of a proliferation of new light states (from objects wrapping the cycle that 
become massless when it collapses). In field theory, the procedure of integrating out these 
states is no longer valid when they become massless, so we should expect that approaching 
this point in the moduli space leads to divergences. In the opposite limit, Vol(S) — > 00, KK 
modes become light, and a similar divergence occurs. In our examples, this will mean that 
moduli have logarithmic Kahler potentials. 

Further explanation of the origin of moduli fields in string compactifications may be 
found in the textbook [60]. An interesting discussion of generic properties of string theory 
moduli spaces that are conjectured to be a requirement of all consistent theories of quantum 
gravity is found in Ref. [61]. 
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